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Abstract. Motivated by the recent progress in engineering artificial non- 
Abelian gauge fields for ultracold fermions in optical lattices, we investigate 
the time-reversal-invariant Hofstadter-Hubbard model. We include an additional 
staggered lattice potential and an artificial Rashba-type spin-orbit coupling 
term available in experiment. Without interactions, the system can be either 
a (semi)-metal, a normal or a topological insulator, and we present the non- 
Abelian generalization of the Hofstadter butterfly. Using a combination of 
real-space dynamical mean- field theory (RDMFT), analytical arguments, and 
Monte-Carlo simulations we study the effect of strong on-site interactions. We 
determine the interacting phase diagram, and discuss a scenario of an interaction- 
induced transition from normal to topological insulator. At half-filling and 
large interactions, the system is described by a quantum spin Hamiltonian, 
which exhibits exotic magnetic order due to the interplay of Rashba-type spin- 
orbit coupling and the artificial time-reversal-invariant magnetic field term. We 
determine the magnetic phase diagram: both for the itinerant model using 
RDMFT and for the corresponding spin model in the classical limit using Monte- 
Carlo simulations. 
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1. Introduction and time-reversal-invariant Hofstadter-Hubbard model 

Experiments with ultracold atoms in optical lattices can provide novel insight into 
traditional strongly correlated condensed matter phases in real materials [1, 2]. 
Examples are the Mott-insulator to superfluid quantum phase transition in both 
the bosonic [3, 4, 5] and the fermionic Hubbard model [6, 7], the physics of 
the BEC-BCS crossover between a Bose-Einstein condensate (BEC) of bosonic 
molecules and the Bardeen-Cooper-Schrieffer (BCS) pairing state of its fermionic 
constituents [8, 9, 10, 11, 12], or the emergence of quantum magnetic order in itinerant 
systems for large on-site interactions [13, 14, 15]. Cold-atom lattice setups offer various 
advantages: they are free from disorder, can be experimentally probed with single-site 
precision (in two dimensions at least) [5, 16], and they allow unique control over its 
system parameters such as lattice geometry or the nature of the constituent particles. 
Particularly attractive is the precise tunability of the interaction strength between 
the particles via Feshbach resonances, which allows for a quantitative study of the 
role of strong correlations in condensed matter systems. One central goal in the field 
is thus the investigation of possible high-temperature (i-wave superfluidity in the 2D 
Fermi-Hubbard model [17, 18]. 

For a long time, however, the toolbox for quantum simulations was missing 
one important aspect: the equivalent of orbital magnetism. Since cold- atoms are 
charge neutral, they do not experience the Lorentz force due to external magnetic or 
electric fields. One way to simulate magnetism for cold atoms is to set the complete 
system in uniform rotation around a fixed axis [19, 20]. This approach, however, is 
severely limited in experiment as it requires, for example, axial symmetry which is not 
present in optical lattice experiments. Recently, it was experimentally demonstrated 
that artificial gauge potentials for cold atoms may also be generated using laser 
beams [21, 22, 23]. The atom experiences an effective gauge field due to light-induced 
couplings between different atomic internal states. This leads to a dressed basis of 
states, and if the atom-light coupling is spatially modulated the moving atom picks up 
a geometric phase [24, 25]. Importantly, this scheme is applicable in the presence of 
optical lattice potentials, where the hopping element acquires a phase factor from the 
well-known Peierls substitution [26, 27, 28, 29, 30, 31]. Furthermore, these artificial 
gauge fields may easily be generalized to the non-Abelian case or to represent spin-orbit 
interaction by coupling internal atomic levels that correspond to different pseudo-spin 
states of the atom [32, 33, 34, 35, 36, 37, 38, 39, 40, 41]. The effective electromagnetic 
fields induced by these artificial gauge fields can be large, and in certain cases they can 
even give rise to the quantum (spin) Hall effect [42, 43, 44, 45], thus allowing for its 
study in a completely new experimental setting [46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56]. 

Here, we investigate the role of interactions in a spinful and time-reversal-invariant 
generalization of the paradigmatic Hofstadter model [57] using real-space dynamical 
mean-field theory (RDMFT) and classical Monte-Carlo simulations [58, 59, 60, 61]. 
We complement our numerical results with analytical arguments. This time-reversal- 
invariant Hofstadter-Hubbard model [62] can be realized using cold- atoms in artificial 
gauge fields [51, 63]. It hosts a number of fundamental physical problems that we 
address below such as the stability of topological phases with respect to interactions 
and the development of exotic magnetic order at large interactions due to spin-orbit 
coupling. 
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The original Hofstadter model considers spinless fermions confined to a two- 
dimensional square lattice potential and subject to a magnetic field in z direction 
perpendicular to the plane. It describes electrons in a magnetic field giving rise to 
the (integer) Quantum Hall Effect (QHE). The single-particle spectrum as a function 
of the magnetic field strength was first computed by D. Hofstadter and exhibits the 
shape of a butterfiy [57]. For a rational magnetic fiux per plaquette a = p/q, in 
units of the Dirac fiux quantum = h/e and with coprime integers the system 
remains translationally invariant with an enlarged unit cell of q lattice sites. All bulk 
bands carry non-zero Chern number C^, and in all gaps one finds chiral edge modes. 
By the bulk-boundary correspondence, their total number is given by the sum of the 
Chern numbers of the filled bands below the gap: Nc = ^^n^^ Ci [64, 65]. For even 
values of q the middle bulk bands touch each other q times at zero energy, where they 
form linearly dispersing Dirac cones [66, 67]. 

In contrast, here we consider the time-reversal- invariant Hofstadter-Hubbard 
model of spinful and interacting fermions. To restore time-reversal symmetry we 
imagine applying a magnetic field in the z direction that only couples to the |t) -spins. 
A second field of the same strength that only couples to the ||)-spins is applied in 
the opposite {—z) direction. We thus obtain a spinful and time- reversal- invariant 
version of the fundamental Hofstadter problem. While such a scenario is certainly 
not feasible with a "real" magnetic field, remarkably it can be experimentally realized 
using artificial gauge fields in cold atoms. The semi- metal (SM) at even values of q 
now becomes a generalization of graphene with q Dirac cones. The number of Dirac 
cones can be tuned by the artificial magnetic field strength and may become large. 
Moreover, the energy gaps which were crossed by Nc chiral edge modes per edge in 
the QHE setup are now traversed by = Nc Kramer's pairs of helical edge modes 
per edge. The case of odd N^ thus corresponds to the topological Quantum Spin Hall 
(QSH) phase, while an even number of A^^ corresponds to a normal insulator (NI) 
phase [65, 43]. 

Particles with opposite spin may occupy the same lattice site and interact with 
each other. The interaction in cold-atom systems is short-ranged. We thus consider a 
Hubbard on-site interaction term in the Hamiltonian. The interaction strength may 
be tuned by Feshbach resonances or by adjusting the lattice depth. Apart from the 
time-reversal magnetic field and the Hubbard interaction term, there are additional 
terms available in the cold- atom experiment [51, 63], such as a staggered optical 
lattice potential and a Rashba-like spin-orbit coupling term that breaks the axial 
spin symmetry as it induces spin fiips between states |t) and |^) when the particle 
tunnels from site to site. 

The Hamiltonian of the time-reversal- invariant Hofstadter-Hubbard model with 
the additional terms takes the form 

j 

(1) 

Here, = (cj^^Cj^) is a spinor of the fermionic creation operators at lattice site 
j = (x, y) with integer x, ;r = (1, 0), y = (0, 1) are unit vectors of the square lattice 
and <j^'^ are the usual Pauli matrices. The parameters tx {ty) denote the hopping 
amplitude in the x-{y-) direction. We focus on isotropic hopping tx = ty = t here, 
and express all energies in units oft = 1. The value of a determines the strength of the 



Correlated Topological Phases and Exotic Magnetism with Ultracold Fermions 4 

(artificial) magnetic field for either spin species which penetrates a lattice plaquette 
in units of the Dirac flux quantum. The on-site interaction strength is denoted by [/, 
the parameter A^^ > introduces a staggering of the optical lattice potential along the 
x-direction, and the continuous parameter 7 controls the relative amplitude of normal 
and spin-flip tunneling in the x-direction. Non-zero 7 breaks the axial spin symmetry 
and introduces mixing between the spin states , Since [(7^,(7^] ^ the 

applied gauge field configuration is non-Abelian if both a,7 > 0. For 7 = n/2 with 
n G Z the part that is off-diagonal in spin space, which is proportional to sin27r7, 
vanishes. We may thus restrict 7 to the interval [0,1/2]. For U = the model in 
Eq. (1) was studied in Ref. [51]. 

The remainder of the paper is structured as follows: in Sec. 2, we first study the 
energy spectrum and the phase diagram of the non-interacting system as a function of 
the parameters a, 7, A^^. We present the equivalent of the Hofstadter butterfly for non- 
zero spin- mixing 7 > and staggering A^^ > 0. Then in Sec. 3, we investigate the role 
of the Hubbard interaction in detail using RDMFT. We present the interacting phase 
diagram in Sec. 3.1, and explicitly determine the robustness of the topological phase 
with respect to the Hubbard-t/. In Sec. 3.2, we show that interactions may drive 
a transition from a trivial to a topological phase. Finally in Sec. 3.3, we discuss 
the magnetic order exhibited by the system at large interactions and half filling. 
The magnetic phase diagram we observe is rich, containing various quantum phase 
transitions due to the TR-invariant magnetic field term a and the spin-mixing term 
7. Both terms may be tuned experimentally. Apart from Neel and collinear magnetic 
order we find various spiral phases. To understand the magnetic phase diagram we 
analytically derive the effective quantum Heisenberg Hamiltonian in this regime. We 
analyze this Hamiltonian in the classical limit using Monte-Carlo simulations and 
recover the same phases that we found in our RDMFT analysis of the itinerant model 
at large U. We end with a conclusion in Sec. 4. 

2. Non-interacting phase diagram: spectrum and non-Abelian Hofstadter 
butterfly 

In this section we consider the non-interacting system and set U = 0. For 7 = 
0, the two spin species decouple and we obtain two time-reversed copies of the 
Hofstadter problem [57]. Fermions acquire a magnetic phase that depends on x 
when tunneling along the ^/-direction. Opposite spins acquire opposite phases due 
to the factor of a^. The system remains periodic on the lattice if we choose the 
magnetic flux per plaquette to he a = p/q with coprime integers p^q. The unit cell 
contains q sites, and consequently the magnetic Brillouin zone (MBZ) extends between 
{kxiky) e [— ^, ^) X [— 7r,7r). The energy spectrum consists of q bands. The bands 
are separated by finite energy gaps, except for the two middle bands in case of even 
(7, which touch each other at q distinct wave- vectors in the MBZ, where they form 
Dirac cones. Specifically, this occurs at /c = {—^^ ±(2r^+i)zr^ for = 2 + 4n and at 
k = (0, for (7 = 4 + 4n with n, m G {0, 1, . . .}. Considering periodic boundary 

conditions (PBC) along the ^/-direction and open boundary conditions (OBC) along 
X, i.e., for cylindrical geometry, we present the energy spectrum for a = 1/10 in 
Fig. 1(a). Within the gaps of the bulk spectrum (blue), we find a varying number of 
helical edge modes (red), specifically = {1,2,3,4,4,3,2,1} from bottom to top. 
This can be understood from an analogy to the corresponding QHE problem where 
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Figure 1. Energy spectrum for a cylinder geometry as a function of ky. Bulk 
spectrum is shown in blue, edge states in red. Artificial magnetic fiux is o; = 1/10, 
and 7 and Ax take values as indicated. Doubly degenerate helical edge states (red) 
are traversing the bulk gaps. 



the time-reversal invariant magnetic field is replaced by a real magnetic field, which 
breaks time-reversal symmetry. There, each bulk band carries a Chern number of 
C = 1, except the two touching bands around zero energy which in total carry a 
Chern number of C = —8. Indeed, the number of chiral edge modes per edge is 
given by the sum over the Chern numbers of the filled bands Nq = fiHed 

d. The 

sign of Nc denotes the direction of the chiral edge mode. In presence of time-reversal 
symmetry, the sum over the Chern numbers for each gap of course vanishes, since up 
and down spins move in opposite directions and thus ^^led ^1 ~ ~ Xli filled ^1 • 

The bulk spectrum in Fig. 1(a) is invariant under translation ky ^ ky = ky-\-27r/q. 
This follows from the fact that for 7 = A^^ = the Hamiltonian H in Eq. (1) commutes 
with the magnetic translation operator which is defined by T^cj = cj"^^ cj-^xTx 
with (J = exp(i27ra) [68, 69]. The phase factor cj"^^ accounts for the gauge field that 
breaks translational symmetry along x. The magnetic translation operator Ty remains 
TyCj = Cj-^yTy. Notc that while [H^Tx] = [H^Ty] = 0, if 7 = A^^ = 0, the translation 
operators do not commute between themselves [T^^Ty] ^ 0, but only [T^^Ty] = 0. 
From TxCk^^ky = Ck^,ky+2TTa follows the additional degeneracy of the bulk spectrum 
(see Fig. 1(a)). 

Let us now investigate the effect of non-zero 7, A^^ 7^ on the spectrum. If either 
of the two is non-zero, the Hamiltonian does not commute with Tx anymore. While 
= ^j{—'^Y^xc\cj still commutes with even powers of Tx^ i.e., [Hx^^T^'^] = for 

n G N, one finds that the ta^-hopping term only commutes with Tx^'^ for < 7 < 0.5. 
This follows directly from 

[^-' = 4+^^^" sin(2^7) [cos(4^a^) - 1 - ia' sm{47ray) 

3 3 

(2) 

As an example, we present results for a = 1/10 with non-zero 7 = 0.25 and A^^ = 
in Fig. 1(b). We observe that the bulk spectrum is still invariant under the shift 
ky ^ ky ^ 27ra| = ky -\- ir. In contrast, if both 7, A^^ 7^ as in Fig. 1(c), the bulk 
spectrum has lost any translational symmetry. Note that it is still doubly degenerate 
at each /c-point and symmetric under {kx^ky) {—kx^—ky) which follows from time- 
reversal invariance and inversion symmetry. 
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Figure 2. Hofstadter butterfly: energy spectrum (black) as a function of 
magnetic flux a = p/q. We sample 500 different rational flux values a = 
{1/2, 1/3, 14/46, 15/46}. In (b) and (c) we observe that the system is 
predominantly metallic. The stripy appearance is due to the fact that we take a 
to be a rational number. The color inside the gap denotes QSH phase (blue) and 
NX phase (green). Panel (a) is for 7 = A^c = 0, (b) for 7 = 0.1, Ax = 0.5, and (c) 
is for 7 = 0.25 and A^ = 1.0. Note the QSH phase (blue) that appears at Ep = 
for q; = 1/6 in panel (c). 



The non-interacting system can thus be either (semi-) met ahic, if the Fermi energy 
Ep hes within a band or at a Dirac point, or it can be insulating, if Ep hes within a 
gap. The normal and topological insulator are distinguished by a Z2 index v which 
is equal io v = mod 2, where is the number of helical Kramer's pairs per 
edge. While we compute the Z2 index by counting the number of edge states for the 
interacting system, we use the efficient method introduced by Fukui and Hatsugai [70] 
in case oiU = 0. 

With the method of Ref. [70] of calculating the Z2 invariant, we have obtained 
the non-Abelian generalization of the Hofstadter butterfly for various values of 7 and 
\x. The results are shown in Fig. 2. In Fig. 2(a), we show the familiar Hofstadter 
butterfly for 7 = A^^ = 0, which shows the spectrum, i.e., the bands, (in black) as a 
function of flux a = p/q. Within the gaps, we have computed the Z2 index and we 
denote topological (normal) insulator by blue (green). In Fig. 2(b), we present the 
non-Abelian butterfly for 7 = 0.1 and A^^ = 0.5. Here, the translational degeneracy 
of the spectrum is lifted and the system is metallic in most parts of the energy range 
(see also Fig. l(b-c)). Some gaps with v = 1 have remained, however, for a < 1/2. In 
Fig. 2(d), we show the butterfly for maximal spin- mixing 7 = 0.25 and A^^ = 1. New 
gaps have opened, and in particular, we now flnd a TI phase at half-fllling Ep = ^ 
for a = 1/6. Similarly, we flnd u = 1 for a = 1/10 and slightly larger (and smaller) 
values of A^^. To demonstrate, we choose two values of a = 1/6 and a = 1/10 and 
show the full Xx-J phase diagram at half-fllling in Fig. 3. The QSH phase appears 
around maximally spin mixing. Since the gap size is larger for a = 1/6 this value 
is experimentally preferable. The case of half-fllling is particularly interesting, since 
interactions drive a transition from a topological to an insulating state with exotic 
magnetic order, as we discuss in Sec. 3. 

Although we have found the QSH phase at half-fllling only in the presence of 
spin-mixing 7 > 0, it is not a requirement to break axial spin symmetry [71, 72]. 
Another way of opening the gap and inducing a QSH phase at half-fllling without 
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Figure 3. Non-interacting phase diagram for a = 1/6 and o; = 1/10 at half- filling 
Ep = as a function of 7 and A^. The QSH effect phase appears symmetrically 
around maximally spin-mixing 7 = 0.25. For our finite size calculations, we 
identify the metallic regions when the "gap" size is smaller than the level spacing. 
Region of QSH effect and gap sizes are larger for a = 1/6 than for a = 1/10. 

breaking axial spin symmetry is to add a diagonal hopping term of the form 

Hd = -tdE,(sW^e''"""'^"^'^'^+S---y^''"""'^"^ [^^]- Our choice of 

Rashba-type hopping 7 and staggering A^^ is thus only an experimentally motivated 
way [51] of finding a QSH phase at half-filling. 



3. Interaction effects: correlated topological phases and exotic magnetism 

One of the key advantages of simulating complex Hamiltonians with cold atoms in 
optical lattices is the ability to control the strength of the interaction via Feshbach 
resonances. Strong interactions and topologically non-trivial band structures can both 
give rise to novel phases of matter. In particular, when interactions and topological 
band structures are simultaneously present, unexpected exotic phases might arise, 
and novel phenomena emerge due to the interplay of topology and electron-electron 
interactions. Recently, the effect of on-site interactions on topological band structures 
has been extensively studied [74, 75, 76, 77, 78, 79, 80, 81]; for details we refer the 
reader to a recent review article [82]. 

Here, we shed light on this issue using RDMFT and thoroughly study the effect of 
interactions on the band structure discussed above. We explicitly show the robustness 
of the edge states up to relatively strong interactions. Further, we identify a scenario 
where the interactions drive the system from a normal to a topological phase. At half- 
filling, the system exhibits magnetic order at large /7, and we find the magnetic phase 
diagram using RDMFT. We also rigorously derive a quantum Heisenberg Hamiltonian 
in this regime, and investigate it numerically in the classical limit of large spin S. 



3.1. Interacting phase diagram and stability of topological phases 

Let us now study the influence of on-site interaction on the phase diagram. We 
use RDMFT to determine the spectral function A{x.,ky^uj) in a cylinder geometry 
with PBC along y and OBC along x. This allows us to directly map out and count 
the helical edge states, which determines the Z2 index v. Edge states are also used 
experimentally to measure v [51, 83, 84, 85]. In the following we focus on fixed 
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a = 1/6, where all relevant phenomena that occur in this system for general p/q are 
qualitatively captured. 

Away from half-filling, we have investigated the stability of the topological phases. 
Placing the Fermi energy in the lowest gap at filling factor np = 1/3 fermions per 
site, we find that within RDMFT these phases appear to be robust with respect to 
interactions. We have computed the spectral function for a finite size system using 
RDMFT and the helical edge mode is clearly visible even for U = 10. 

At half-filling, we have determined the interacting phase diagram for a = 1/6 as 
a function of A^^ and U for various values of 7, and it is shown in Fig. 4. Since the 
phase diagram is symmetric around 7 = 0.25 is it sufficient to consider < 7 < 0.25. 
At U = the system is a (semi)-metal for 7 < 0.2 and A^^ < 1.26. It turns into a 
normal band insulator for A^^ > 1.26 (see Fig. 3 (left)). At large U > Udj^Xx) the 
system exhibits magnetic order. The type of magnetic order depends on 7, as we 
show in detail in Sec. 3.3. The critical U increases with A^^, which follows from the 
fact that Xx tends to localize fermions in every other column whereas U favors a state 
with exactly one fermion per site. In general, increasing the interaction strength U 
reverses the effect of the staggering potential A^^. This leads to the interesting situation 
at 0.2 < 7 < 0.3, where we find interaction driven transitions from a metallic or a 
normal insulating phase into the topological QSH phase. The QSH phase then extends 
over a broad range of about U^2 — A before the system eventually develops magnetic 
order at /7 = 4 — 5 and becomes gapped again. 
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Figure 4. Interacting phase diagram at o; = 1/6 as a function of interaction U 
and staggering for various values of spin-mixing 7. We find (semi)-metallic 
regime (Metal), magnetically ordered state (Mag), normal insulator (NI) and 
quantum spin Hall state (QSH). 



3.2. Interaction induced topological quantum phase transition 

The non-trivial topology of topological insulators is an effect entirely due to the band 
structure, i.e., the topology is already contained in the non-interacting part of the 
Hamiltonian. On the other hand, it was realized that topological band properties 
can also arise in topologically trivial systems from the spontaneous breaking of a 
symmetry due to interactions [86] . This idea suggests that topological insulators could 
be present in a much larger class of materials. Recently a few theoretical models have 
been proposed where this scenario is realized [62, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95]; 
from the experimental point-of-view, the situation is, however, less satisfying and an 
experimental study of an interaction-driven topological insulator is still lacking. A 
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Figure 5. Intensity plot of the single-particle spectral function A(ky,uj) = 
fj^ Aikx^ky^uj) as a function of uj and ky for a = 1/6, = 1.5, 7 = 0.25 
and different interaction strength U . Left panel is for U = 0.5, where system is 
in the normal insulating phase, (middle) is for U = 1.0 where we find a metallic 
phase, and (right) is for U = 3.0, where the system is in the quantum spin Hall 
phase. The spectral function visualizes the topological difference, in (left) no 
edge state is crossing the bulk gap while in (right) a pair of helical edge states is 
traversing the bulk gap. 

simple example of an interaction-induced topological quantum phase transition into 
a QSH phase is realized in the Hamiltonian (1) for a = 1/6, i.e., here we provide an 
example which is experimentally relevant. We will discuss this topological quantum 
phase transition in the following. 

In principle, interactions are encoded in the self energy S(/c,cj). In most 
approaches, one needs to approximate the self energy for practical reasons. Usually 
one uses approximations where the self energy becomes either purely momentum- 
dependent or purely frequency-dependent. 

In static mean-field, for instance, the self energy is reduced to (or even to 

a constant Eq) and enters as a band structure renormalization. In certain cases, 
such a band structure renormalization can cause band-inversion and induces the 
topological phase. Another proposal discussed a fiuctuation-induced topological phase 
transition [88, 95], i.e., a topological phase induced solely by T^{uj) being momentum- 
independent. 

In our Hamiltonian (1), the topological phase transition can be understood as a 
So-driven transition. This constant part of the self energy is contained as a lowest 
order contribution in both and I](cj). The latter quantity is available in DMFT- 
like methods and thus we can observe this interaction-induced topological phase 
transition. The transition appears for parameters such as 7 = 0.22 and A^^ = 1.5 
where the band structure is in a topologically trivial phase. The reader should note 
that large values of the staggered potential A^; avoid the QSH phase, i.e., for the same 
parameters but A^^ = 0.75 the QSH phase would be present. By increasing U we 
effectively reduce the effect of A^^ and drive the system into a metallic phase. Under 
further increase of U the gap re-opens and one enters the QSH phase. Consequently, 
these transitions from trivial insulator to metal and from metal to QSH phases are 
due to interactions which cause a band structure renormalization. We have shown the 
corresponding spectral functions in Fig. 5 for the three different regimes. Absence or 
presence of edge states, respectively, distinguish the trivial and topological insulator 
phases. 

3.3. Tunable exotic magnetism 

We learned previously that at half-filling, the system turns into a magnetically ordered 
insulator at a = 1/6 and large interaction strength U ^ t (see Fig. 4). This behavior 
occurs generally in the half-filled TR-invariant-Hofstadter-Hubbard model at large 
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Figure 6. Magnetic order of the itinerant model at t/ = 5 as a function of 
7 G [0, 0.25] for different values of a = 1/3, 1/6, 1/4 (from top to bottom) obtained 
within RDMFT. We generally find Neel order around 7 = and collinear order 
around 7 = 0.25. In between for 0.05 < 7 < 0.20, we observe different kinds 
of spiral order with periodicities of three, four, six, and eight sites due to a 
competition of a Dzyaloshinskii-Moriya (DM) and an anisotropic XXZ Neel type 
spin interaction, see Eq. (4). Phase boundaries are indicated by solid vertical 
lines. Black dots denote the actual values of 7 that we simulated using RDMFT. 
The magnetization diagrams show the direction of magnetic order in Sy-Sz- 
plane. The background color indicates the magnetization amplitude Stot{x). The 
magnetization diagrams are extracted from the specific parameter point that is 
indicated by the dashed circle. If the classical spin Hamiltonian exhibits a different 
phase boundary than the itinerant model, we denote the phase boundary of the 
spin Hamiltonian by a vertical dashed line with label ii^spin- 



U . Exotic magnetic order appears that depends on the value of the Rashba-type 
spin-orbit coupling 7 and the artificial magnetic fiux a. Importantly, by tuning 
these parameters one crosses various magnetic quantum phase transitions from Neel 
order to various spiral ordered phases to collinear order. Similar phases have been 
reported at strong interactions in other cold-atom setups that include non-Abelian 
gauge fields [96, 97, 98, 99]. 

Using RDMFT we determine the magnetic phase diagram in Figs. 6 and 7 as a 
function of 7 for various values of a. We work at fixed interaction strength [7 = 5. 
We complement our analysis by rigorously deriving a quantum spin Hamiltonian at 
small t'l^y/U when charge fluctuations freeze out and we numerically determine its 
magnetic phase diagram in the classical limit of large S using classical Monte-Carlo 
simulations, and compare to the RDMFT results of the itinerant model. While the 
phase boundaries are often slightly shifted, we generally find the same type of magnetic 
order in both the itinerant model at = 5 and the classical spin Hamiltonian. 

In Figs. 6 and 7, we observe that the order along the y direction is always 
antiferromagnetic, and we find various types of order along the x direction. In all 
cases investigated, the Sx component is zero, and we can represent all spins in the 
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Figure 7. Rich magnetic phase diagram at a = 1/5, f/ = 5 as a function of 
7. See figure 6 for notation. We note that due to quantum fluctuations, the 
boundaries between different types of order for 0.05 < 7 < 0.125 in the itinerant 
model, which are indicated by vertical solid lines, is signiflcantly shifted when 
compared to the phase boundaries of the classical spin Hamiltonian, which are 
indicated by vertical dashed lines with text ii^spin- 



Sz — Sy plane. For < 7 < 0.05 one finds usual Neel order for all four values of a 
that we consider. Around maximal spin-mixing 0.2 < 7 < 0.25 the order is collinear 
with spins ordering ferromagnetically along the x-direction and antiferromagnetically 
along the ^/-direction. In between, the system exhibits various forms of spiral order 
with different periods of three, four, five, six, eight or ten lattice sites for the values 
of a that we present here. 

We can describe the spiral states using a variational ansatz with only four 
parameters. The order of the spin Hamiltonian for each set of parameters is uniquely 
defined by the angle of rotation Oi between nearest-neighbor spins along the x- 
direction, i.e., cos{6i) = Si • Si-^x/\Si\\Si-^x\' A single ansatz of the form 

= ^0 + ^u; cos(a;x + r]) (3) 

with the four variational parameters {(/)o, (/>u;, cj, 7^} perfectly fits the order found. Note 
that 6i = 00 is sufficient to describe Neel, collinear, and simple spiral order, which 
are indicated in Figs. 6 and 7 by red, yellow and orange regions respectively. We 
distinguish these from the more complex spiral order, indicated by green regions, that 
can only be described with one additional Fourier component of frequency uj. We 
determine these variational parameters by a Fourier analysis of the Monte-Carlo data 
and, as one might expect, these parameters are close to nice analytical values. For 
example, a = 1/6, 7 = 0.15 has the values 0o ^ 7r/3, u ^ 27r/3, (j)^ ^ 7r/6. 

The magnetic phase diagram is symmetric around the point 7 = 0.25, except that 
the spiral phases have different chirality for 0.3 < 7 < 0.45. Quantum fluctuations 
tend to reduce the amplitude of the magnetization 5'tot < 1/2, which depends on 
the values of 7, a and U. In addition, for collinear and some spiral orders 5'tot(^) 
is also spatially staggered along the x direction for intermediate interaction strength 
[/, which we indicate by the colored background of the spin boxes in Figs. 6 and 7. 
As shown in Fig. 8, this spatial staggering reduces with increasing U as the fermions 
become more and more localized. The type of magnetic order that occurs can be 
understood by rigorously deriving a quantum spin Hamiltonian at large interactions 
when the charge degrees of freedom freeze out and only virtual hopping events take 
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where Ji = ^tl/U and the spin operators are defined in terms of the fermionic 
operators as Sj = ^c^-crcj with cr = (a^, a^, a^). The first part describes spin exchange 
in X direction, which contains the Rashba-couphng 7. For 7 = n/2 with n G Z, the 
couphng takes the form of the famihar Heisenberg Hamiltonian. For other values of 7, 
however, the SU (2) symmetry is broken down to XXZ-type with anisotropy direction 
5*^ in spin space. For 7 7^ n/4 an additional Dzyaloshinskii-Moriya (DM) interaction 
term in the YZ plane is present. It is this term that is responsible for the spiral order 
that we find for intermediate values of 7. In the absence of an artificial magnetic field 
{a = 0) no direction is preferred for spiral order (in YZ) or antiferromagnetic order 
(in X). 

The second part of Eq. (4) describes spin exchange in y direction, and depends 
on the artificial magnetic field strength a. It is periodic along the x direction with an 
extended unit cell of ^^72 lattice sites. For a = 1/2 one obtains the ordinary Heisenberg 
exchange, while for other values of a the spin interaction of the 5*^ and components 
exhibit an x dependent modulation of the variable strength of Heisenberg exchange 
and DM-type interaction. The spin interaction along Z is always antiferromagnetic. 
In the absence of a Rashba-like term (7 = 0) the periodic modulation causes ordering 
in the X and Y directions to be unfavorable, such that only Z antiferromagnetic order 
occurs. For 7 7^ 0, the competition between the present terms causes the X direction 
to be unfavorable and only order in the Y-Z plane is observed. 

To make a connection between the spin model (4) and the itinerant model (1) at 
large [/, we have solved for the phase diagram of the spin model in the classical limit 
of large spin S using Monte-Carlo simulations. As shown in Figs. 6 and 7, we find 
that while the phase boundaries between the different magnetically ordered regions 
are shifted in some cases, which is not surprising considering that we are taking the 
classical limit of a 5* = 1/2 quantum spin model, the type of magnetic order that 
occurs in the two models is the same. 



4. Conclusions 

In conclusion, ultracold atoms in optical lattices in the presence of non-Abelian gauge 
fields provides a rich platform to study the influence of interactions on topological 
phases. Motivated by experiment, we have investigated a specific setup, the time- 
reversal- invariant Hofstadter-Hubbard model, and determined the interacting phase 
diagram. Interactions drive various phase transitions. We identified a scenario where 
a normal insulator becomes topologically non-trivial for increasing interactions. At 
large interactions and half-filling, the systems develops tunable exotic magnetic order 
with Neel, collinear, and spiral phases due to the interplay of Rashba-type spin- 
orbit coupling and the artificial magnetic field. We explicitly prove the robustness 
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Figure 8. Spatial staggering of the magnetization amplitude Stot{x) as a function 
of U for q; = 1/6. Staggering is defined by the difference of the maximal and the 
minimal S'tot within one period of the spin order. The staggering is indicated 
by the colored background in the spin boxes of Figs. 6 and 7. We observe that 
the staggering decreases for increasing U, because the fermions become more and 
more localized and quantum fluctuations are suppressed. This also occurs for the 
other values of a that we investigated. 

of helical edge states in the correlated topological phases, which is crucial for cold- 
atom detection schemes. 
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